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ABSTRACT 

For a number of starless cores, self-absorbed molecular line and column density observations have implied the 
presence of large-amplitude oscillations. We examine the consequences of these oscillations on the evolution of 
the cores and the interpretation of their observations. We find that the pulsation energy helps support the cores and 
that the dissipation of this energy can lead toward instability and star formation. In this picture, the core lifetimes 
are limited by the pulsation decay timescales, dominated by non-linear mode-mode coupling, and on the order of 
~ few x 10 5 -10 6 yr. Notably, this is similar to what is required to explain the relatively low rate of conversion of 
cores into stars. For cores with large-amplitude oscillations, dust continuum observations may appear asymmetric 
or irregular. As a consequence, some of the cores that would be classified as supercritical may be dynamically 
stable when oscillations are taken into account. Thus, our investigation motivates a simple hydrodynamic picture, 
capable of reproducing many of the features of the progenitors of stars without the inclusion of additional physical 
processes, such as large-scale magnetic fields. 

Subject headings: hydrodynamics-ISM: clouds - ISM: globules - stars: formation 



1. INTRODUCTION 

The starless cores owe their reputation as the sites of 
future star formation to the similarities they share with 
neighboring cores already in the pains of stellar birth. Small 
(tenths of pc), dense (« H , ~ 10 3 to 10 6 cm -3 ) dark clouds 
of a few solar masses, the starless cores contain no infrared 
sources above the sensitivity level of the IRAS satellite 
(about 0.1 L Q at the distance of Taurus). Mixed in with a 
population of cores, embedded infrared protostars, and young 
T-Tauri stars , the starless cores seem on t h e verge of issue. 
(IMvers et alj 119831: iMvers & Benson! 119831: iBenson & Mversl 
1989; Beichman et al.l 1 19861 Idi Francesco et al J 120071 : 



Ward-Thompson et al.ll2007t IBergin & Tafallall2007l) 



The starless cores are remarkable in another sense. They may 
be the only structures in the interstellar medium in quasi-static 
equili brium with Ufetimes l onger than their sound-crossing 
times (|Beichman et al.lll986HJessop & Ward-Thomp son 2000; 
iKirk et alj 120051: IWard-Thompson et all 120071) . Whereas all 
larger-scale clouds are better d escribed as trans ient struc- 
tures within a turbulent cascade dField et alj 120081) . the star- 
less cores show observed dens i ties t h at are cons i stently 
matched (IWard-Thompson et al.l [19941: lAndre etall 119961: 
IWard-Thompson et al.il 19991: iBacmann et alJ2000tlEvans et al. l 
2001 1 : lAlvesetal.l 120011 : iTafaUa et al.l 120021: iKandori et alJ 
2005: lKirk et al.l2005l) by the truncated solutions of the isother- 



mal Lane-Emd en equation, i.e., the Bonnor-Ebert (BE) spheres 
(Bonnor 1956). These are self-gravitating spheres supported 
by thermal energy and bounded by an external pressure. The 
source of the external pressure appears to depend upon the 
star forming region, though its necessity is well documented. 
Within the Pipe nebula, the core energetics are consistent with 
a constan t external pressu re supplied by hotter, more rarei- 
fied gas dLadaet al.l 120081) . In other regions, such as Lu- 
pus (Teixeira et al.l 120051) . the cores appear to be surrounded 
by lower density molecular gas with observed molecular line 
widths > 1 kms. Here and in similar regions small-scale 
(smaller than the observing beam) supersonic turbulence in the 
low-density molecular gas around the cores may supply the 



confining pressure. 

While a continuous extension to the larger-scale ISM 
prompts the question of whether the cores are individual dy- 
namical entities, a distinct boundary is evident in several ob- 
servat ions. First, maps of dust absorption in the Taurus re- 
gion (IBacmann et al.l 12000) often show sharp edges defin- 
ing the core boundaries. Second, around at least one core, 
B50, a boundary is defined by a sharp transition (within less 
than an observing beam width) between subsonic and super - 
sonic turbulence in the molecular gas dPineda et alj 120101) . 
Third, a boundary is a neces sary feature in models of radia- 
tive c o oling and UV heatin g fevans et al.l 1200 It IShirlev et ail 



20021: IZucconi et al.l [200lt IStamatellos & Whitwort 



lev et alj 
hi 120031: 



Goncalves et al.l 120041 iKeto & Field! I200l " that use the ba- 
sic BE structure an d accurately predict the observed gas and 
dust t emper a tures dWard-Thompson et al.l 120021 : iPagani et al.l 
120031 120041: iCrapsi et al.l 120071) , the observed variation in 
the excitation of the CO molecule across the starless core s 
(IBergin et al.ll2006t iSchnee et al.ll2007t IKeto & Casellll2008l) . 
and the observed bri ghtness and profiles of molecular lines 
( IKeto & CaseUill20Toh . 

The predictive success of models based on BE spheres (in- 
cluding a core boundary) implies that the starless cores are 
indeed individual dynamic entities in quasi-static equilibrium, 
balancing the forces of their own self-gravity, internal energy, 
and the larger-scale ISM. In some cases, this equilibrium may 
be unstable, leading the core to contract, though their internal 
pressure is sufficiently high to prevent free-fall. 

The hydrodynamic equations that describe the equilibrium 
configuration of both stars and the starless cores allow for per- 
turbations in the form of sonic waves that may persist for many 
crossing times. Such oscillations are observed in the Sun, other 
main sequence stars, and white dwarfs, and given the turbu- 
lent state of the larger-scale IS M, should be exp ected in the 
starless cores as well. Analytic (IKeto et al.ll2006l) and numer- 

1 The core B5 contains a protostar but is otherwise similar to the starless 
cores. Since the effects of the protostar are not significant at the core boundary, 
the same demarcation would be expected in the truly starless cores. 
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ical dKeto & Fieldll2005t iBroderick et~alll2007l) models of os- 
cillating cores easily reproduce the characteristic asymmetric 
shapes of molecular lines seen in the starless cores. While 
the spectral-line profiles of the low-order modes (the dipolar 
and breathing modes are observational indistinguishable from 
simple contraction, expansion and rotation) a few cores show 
simultaneous inward and outward mo tions (lAguti et alj |2007: 
lLada et aT]|2003b iRedman et al.l l2006) that are consistent only 
with oscillatory modes. Oscillations produce density pertur- 
bations as well, and the dipole and quad rupole modes nat urally 
reproduce the ellipsoidal morphologies dKeto et al.ll2006l) char- 
acteristically seen in observations of the cores (iBacmann et alJ 
l2000t I Ward-Thompson et a!1[l^lKrrFeTaIll2005l) . 

Observations suggest that these oscillations, rather than a 
minor perturbation on quasi-static equilibrium, provide a sub- 
stantial contribution to the dynamical stability and therefore 
are important in the evolution of the cores toward star for- 
mation. First, a significant fraction of starless cores have 
observed density profiles that can only be matched by a BE 
sphere if the gas temperature in the model sphere is made 
considerab ly higher than is known to be the c a se in the star- 
less c ores (lEvans et alJ 1200 It iKirk et all 120051 : iKandori et alJ 
120051) . Second, the sizes of many of the starless cores im- 
ply that they sho u ld be unstable to g ravitational collapse 
(lEvans etaflboOU IKandori et ail 12005b . However, statisti- 
cal analyses of the relative numbers of stars and cores in a 
sample of different star-forming regions imp ly lifetimes that 
are several times the relevant free-fal l time dBeichman et al. | 
1986 1: Uessop & Ward-Thompsonll2000t IWard-Thompson et alJ 
20071) . Reconciling these observations require a source of en- 
ergy within the cores in addition to the thermal energy. In- 
ternal oscillations or turbulence are suggested by the obser- 
vation that the spectral line widt hs in many cores are larger 
than the expected thermal width (|Dickman & Clemens! 1 19831 : 
iMyers & Benson! 119831: lLada et alJ 120081) implying internal 
subsonic velocities on spatial scales unresolved in the observ- 
ing beam. 

These observations combined with modelling suggest the 
following evolution for the starless cores: The cores are formed 
at the bottom of the super sonic turbulent ca scade that domi- 
nates the larger-scale ISM dField et alJl2008l) . producing cores 
with a spectrum of masses and internal turbulent energies. 
Those cores massive enough to be gravitationally unstable im- 
mediately begin contraction to form a star, while those with 
appropriate mass and internal energy become self-supported 
and self-gravitating. These cores are decoupled from the cas- 
cade in the sense that they exist as indivi dual dynamical enti - 
ties independent of the larger-scale flows dKeto & Fieldll2005l) . 
The internal turbulence inherited from the parent flow, and 
present when a core forms, is part of its initial dynamical sta- 
bility but is expected to dissipate over time. Thought of as a 
spectrum of individual modes of oscillation, the turbulence de- 
cays by non-linear mode-mod e coupling (collisions between 
modes) ( IBrode rick et al] 120071) and by trans mission through 
the co re boundary into the larger-scale ISM ( IBroderick et alJ 
2008). The gradual loss of this internal turbulent energy leads 
eventually (on time scales of several crossing times) to gravita- 
tional instability and star formation. 

In this paper we explore the latter half of this evolution. We 
develop a model for the internal turbulence of the cores and 
examine its consequences for the evolution toward star for- 
mation and for the interpretation of dust continuum observa- 
tions. Using the fact that the turbulence can be fully described 
as a collection of oscillating modes, we consider three mod- 
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FIG. 1 . — Surface pressure as a function of radius for a pulsating fixed-mass 
isothermal gas sphere. For reference, the static configuration is shown by the 
dotted line, and becomes unstable at a center/surface density ratio, f , of 14.04. 
The red, blue and green lines shown the surface pressure resulting from oscil- 
lations containing 30% of the unperturbed configurations gravitational binding 
energy for the Kolmogoroy (Equation [7), Flat (Equation |5J and fundamental- 
Kolmogorov (Equation 1 10) energy spectra, respectively. Along the top axis f 
is shown until the critical value is reached (above which the mapping from f 
to R is no longer single valued, as evidenced by the spiral structure in P vs. 
R for f > fait). Along the right axis the corresponding mass at fixed surface 
pressure is shown. Finally, within the shaded region at the top no stable static 
configurations exist. 



els for the turbulent energy distribution, a Kolmogorov spec- 
trum, a flat spectrum, and a spectrum dominated by long wave 
lengths. In Section [2] we compute the effective pressure aris- 
ing from this turbulence. We find that the breathing modes 
are destabilizing, but modes with shorter wave lengths gener- 
ate a supporting pressure. In Section [3] we generate random 
realizations of turbulent cores and find that a typical stable, os- 
cillating core would appear to be dynamically unstable, if the 
observations were fit to a static BE model. In Section |H we 
follow the evolution of the turbulent spectrum as it decays. We 
find that the spectrum of the turbulence evolves toward domi- 
nance by the longer wave lengths. In conclusion, we find that 
we can describe the evolution of the cores toward star forma- 
tion as a purely hydrodynamic phenomenon without the neces- 
sity of support by large-scale magnetic fields (though a tangled 
magnetic field may co-exist in equipartition with the turbulence 
without substantially changing the findings). 

2. THE STABILITY OF PULSATING CORES 

In this section we compute the dynamical consequences of 
large-scale oscillations for the underlying starless core. We find 
it is possible for pulsations both to support the core and to ini- 
tiate collapse, the distinction lying in the relative importance 
of the low and higher-order oscillation modes. It is possible to 
do this with a macrosco pic descriptio n similar to the stability 
analysis of a BE sphere dBonnorl 19561) . requiring only the pul- 
sation frequencies (which depend upon the core structure) and 
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the distribution of energy amongst them. Here we briefly re- 
view how stability may be determined by the macroscopic core 
properties, make explicit the notion of "slow" compression, use 
this definition to compute the response of the pulsations, and fi- 
nally assess the consequences for a number of different mode 
compliments. 

2.1. The Stability of Stationary Cores 

The starless cores in the BE model are held together by a 
combination of their own self-gravity and the external pres- 
sure, and these two inward forces are resisted by the outward 
thermal pressure. The dotted line in Figure Q] shows the equi- 
librium configurations as a function of the external pressure,/ 5 , 
and the core size, R. At large radii, the portion of the curve to 
the right of R/ R cl ,, = 1 , the gravitational force at the boundary 
is weak compared to the confining pressure. If we compress a 
core on this portion of the curve, we find that the internal pres- 
sure increases with the increasing density, and a larger external 
pressure is required to confine the core, dP/dR < 0. If the 
compression is removed, the core will re-exand to its original 
equilibrium. However, if the radius becomes sufficiently small, 
R/Rcrit < 1. gravitational force dominate at the boundary, and 
a smaller external pressure would be required to maintain equi- 
librium as the core is compressed, i.e., dP/dR> 0. In this case, 
the core can not resist further compression and is therefore un- 
stable to contraction. For static isothermal gas spheres, this is 
greatly simplified by the fact that the equilibrium configura- 
tions form a one-parameter family, which may be chosen such 
that the relevant parameter is the ratio of the surface and central 
densities (or, equivalently, pressures), £ dBonnorll956l) . In this 
case the maximum ratio is = 14.04. 

The thought experiment we describe above to define the sta- 
bility of cores employs a sequence of cores at fixed mass and 
varying surface pressure. However, it is perhaps more physi- 
cally appropriate to hold the surface pressure fixed, determined 
by the properties of the surrounding molecular gas, and vary 
the core mass instead. That is, at a given surface pressure, core 
radius and core temperature identify the unique, stable isother- 
mal core mass, M(R) (if an equilibrium configuration exists at 
all). The mass and pressure are related by the scaling relations 
that connect gas spheres with identical £'s. Explicitly, from 
the equations of hydrostatic structure, at a given ( we have 
M oc Cj/i 3 '/ 2 , and thus fixing c s (i.e., fixing the core temper- 
ature) gives M oc P~ l l 2 . The P(R) shown in FigureQ]are com- 
puted at fixed mass, i.e., M = const. Conversely, if we choose 
P = const, we have 



M(R) 
M 



P(R) 



"I 1/2 



(1) 



Thus, the equilibrium core mass M(R) at fixed surface pressure 
and core temperature is identical, up to a non-linear rescaling, 
to P(R). In Figure 1 the right-hand vertical axis is labeled with 
M(R) instead of P(R). Most importantly, P cr jt corresponds to a 
critical mass, M cr ; t , above which no stable configuration exists. 

2.2. Response of Pulsations to Slow Compression 

The situation becomes more complicated in the presence of 
energetically significant pulsations. This is because we must 
not only adjust the underlying gas configuration, but also the 
oscillation amplitudes. In principle, a rapid compression can 
excite oscillations itself, and we must carefully specify the dy- 
namical nature of the experiment we use to define stability. 



More generally, we must carefully identify how "equilibrium" 
configurations are defined along the entire P(R) curve. In the 
static case a similar consideration is already made when we 
assume that the core mass and temperature are fixed; both are 
natural physical conditions for the starless cores. Since the core 
has a boundary, its mass is constanQ. Radiation provides an ef- 
ficient thermal coupling to the surrounding medium keeping 
the core temperature unvarying on a time scal e much shorter 
than the dynamical time scales of the cores dKeto & Casellil 
l2010h . 

Although radiative equilibrium violates the adiabatic condi- 
tion in the strict thermodynamic sense, if the compression is 
assumed to occur slowly relative to the oscillation periods, we 
may still treat it as adiabatic for the purposes of evolving the os- 
cillation mode amplitudes. In this case, what remains fixed are 
the adiabatic invariant^ of the oscillations. Computing these is 
simplified by the fact that we may treat the pulsation modes as a 
collection of uncoupled, one-dimensional harmonic oscillators, 
with the mode amplitudes defined by the radial and angular 
quantum numbers n, I and m (with their standard definitions), 
A„i m , governed by 



d Anltn 
dt 2 



^rdrrAntm j 



(2) 



for a mode dependent angular frequency , uj„i m 
dGingold & Monaghanl [19801: iRathore et all l2003h . For 
this the adiabatic invariant is well known to be E n i m /uj n i m , 
where E„i m = M„/,„w 2 /m A 2 , m /2 is the mode energy (in which 
M„i m is the mode "mass" and computed from the mode 
structure directly). Thus, we may compute the evolution of the 
energy in a particular mode during an adiabatic compression, 
given u>„i m (R). Specifically, from the adiabatic invariant we 
obtain 

dEfiiffi E n i m dco n i m 



dR 



Unlm dR 



However, the presence of an additional component to the core's 
total energy associated with the oscillations, E osc = J2nimE" lm > 
also implies the existance of an additional component to the 
surface pressure. That is, generally, under adiabatic compres- 
sions the oscillations will produce a dynamically generated sur- 
face pressure 



p = . 

1 osc — 



dF 
" dV 



d \ 

dV^ 



1 



E„i,„ dio„. 



III! 



nlm 



nlni 



AttR 2 ^-f LOnlm dR 
nlm 



(4) 



Thus, the pulsation contribution to the surface pressure is de- 
pendent upon doj n i,„/dR and the oscillation energy spectrum 
(which we shall discuss in more detail below). 

In the dicussion thus far, we have not placed any require- 
ments upon the underlying static configuration which de- 
fine the structure of the pulsations. Therefore, we can gain 
some insight into what we might expect from a simple, one- 
dimensional example, a uniform temperature and density cu- 
bical cavity of side-length 2R. In this case the wavelength 
of each oscillation mode, the standing waves of the cavity, is 

2 We allow that real cores may occasionaly gain mass from interaction with 
the surrounding ISM, but this is not a necessary process for this discussion. 

3 The notion of an adiabatic invariant arises naturally within the context of 
the Hamilton-Jacobi formulation of classical mechanics. A full discussion of 
adiabatic invariants and their relationship to the action-angle vari ables of the 
Hamilton-Jacobi theory m ay be found in §1 1-7 of Goldstein 1 1980) or §50 and 
Landau & Lifshitz 1 1976). Of particular relevance are the adiabatic invariants 
of the simple harmonic oscillator, which is presented as an example in both 
references. 
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FIG. 2. — Normalized pulsation frequencies for a number of low-harmonic 
oscillations of the isothermal gas sphere as a function of sphere radius, R. The 
thick and thin black dotted line shows the BT [ and fi~ 3 / 2 , respectively. While 
n = 0-5 and / = 0-5 are shown (the frequencies are degenerate in m), only those 
that may be easily distinguished from the RT 1 asymptotic behavior are labeled. 
Frequencies with / = 0, 1, 2, 3, 4 and 5 are colored red, orange, yellow, green, 
cyan and blue, respectively. 



trivially proportional to R and corresponds to a frequency pro- 
portional to R~ . Thus we find dui n i m /dR = -R~ 2 and P osc oc 
T? -4 oc V~ 4//3 , i.e., the phonons constitute a T = 4/3 gas. Of 
course, this may be computed directly from the dispersion re- 
lation u> = c s k where c s is fixed. Nevertheless, we may expect 
short-wavelength pulsations to provide an additional compo- 
nent with a nearly identical equation of state to that of small- 
scale, tangled magnetic fields. 

Generally, u„i„,(R) will depend upon the structure of the 
equilibrium core, in particular upon the relative sizes of the 
core and the wave length of the oscillations. For the oscilla- 
tion modes of isothermal cores, we may construct the mode 
structures and compute the mode frequencies directly, shown 
in F igure [2] for a variety of radii (or, e quivalently, £; see, 
e.g., lRathore et al.ll2003l: iKeto et al.ll2006[) . High-order modes, 
i.e., modes with wavelengths much shorter than the core ra- 
dius, don't probe the gas density gradients, and we recover 
u,,i m (R) oc R~ l . In contrast, low-order modes do show substan- 
tial departures from the R~ l dependence, becoming most sig- 
nificant for the fundamental mode (« = 0). For high meridional 
wave number (/ > 1) the fundamental modes have wo/m oc R~ 3 ' 2 , 
i.e., the pulsation frequencies are roughly proportional to the 
Keplerian frequency at the core surface, lok = (GM/R) 3 / 2 , cor- 
responding to an adiabatic index of 3/2. This is not surprising 
as these modes probe the radial structure of the isothermal core. 

The only oscillation mode substantially deviating from either 
the R~ l or R~ 3 ^ 2 behaviors is the breathing mode (n = 0, / = 
and m = 0), responsible for the onset of instability at £ = £ cr it, at 
which which wooo vanishes. Near £ cr jt, duiooo/dR is positive and 
diverges, yielding a formally infinite negative surface pressure. 



2.3. The Oscillation Energy Spectrum 

Once an energy spectrum, E„i m , has been specified, we may 
produce a revised P(V) which includes the contributions from 
the oscillations. This is now an explicit function of the mode 
amplitudes as well as £ cr i t , and thus no longer has a unique 
shape. In setting the E n \ m we will have intrinsically made some 
statement about the mode excitation and evolution. Indeed, as 
we shall see in Section HI the evolution tends to move power 
from small wavelengths to large wavelengths. Despite this, we 
may bound the likely possibilities, and already a number of im- 
portant qualitative statements can be made about the dynamical 
importance of pulsations. 

2.3.1. Large-Amplitude Monopole Oscillations 

The breathing mode is unique in that, in principle, it can pro- 
duce a negative contribution to the surface pressure, diverging 
at £ = £ cr i t . The breathing mode corresponds to a simple expan- 
sion and contraction of the core, similar to changing C in the 
family of equilibrium configurations. On the inward pulse, the 
contraction may result in ( > £ cr ;, and collapse of the core. This 
can potentially destabilize otherwise stable isothermal cores, 
inducing collapse at a lower external pressure (or equivalently, 
smaller core mass). However, the mathematical analysis of the 
consequences of energetic breathing modes specifically, and 
monopole pulsations more generally, are not clear for a number 
of reasons. 

Foremost among these is the failure of the linear descrip- 
tion of the monopole oscillations. That this must happen near 
C = Ccrit f° r the breathing mode is suggested by the fact that 
wooo vanishes at this point, and thus Aooo = ^odo V 2-Eboo / Mx)o 
diverges for finite mode energies (Mooo is not pathological at 
Ccrit)- More importantly, the density perturbations associated 
with the breathing mode diverge, doing so most strongly at the 
core center. Restricting our attention to pulsation amplitudes, 
and thus energies, for which the central density perturbation is 
less than the unperturbed central density (a somewhat permis- 
sive definition of "perturbation") results in fractional changes 
in the supportable external pressure of less than 0.5% at all £, 
and half that at most That is, when we consider mode am- 
plitudes that are indeed perturbative, the dynamical effect of 
the breathing mode is negligible. A similar constraint limits 
the dynamical importance of the higher-order monopole pul- 
sations. Unlike the breathing mode, these generally increase 
the allowed external pressure, as initially anticipated. Never- 
theless, as with the breathing mode, limiting their amplitudes 
to the perturbative regime results in a small correction to the 
stability of isothermal cores. 

What happens in the non-linear regime is not a priori clear, 
and beyond the scope of this paper. However, we might expect 
that cores with oscillations dominated by non-linear monopole 
pulsations are short-lived. If non-linear monopole oscillations 
can indeed be excited in a subset of starless cores, and these 
do induce instability in a fashion similar to that anticipated 
by the linear analysis of the breathing mode given above, then 
these cores will either collapse or disperse on a single dynam- 
ical time. These cannot be the apparently long-lived cores that 
are observed. That is, the observed cores are necessarily those 
without large-scale breathing modes. 

Alternatively, large-amplitude breathing modes may simply 
be uncommon in starless cores. If the pulsations are driven by 
the turbulence from the parent molecular cloud, it may be diffi- 
cult to excite monopole oscillations. The intra-core motions are 
typically subsonic, and at the core scale, the turbulence is there- 
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fore nearly incompressible. However, the monopole modes 
generally, and the breathing mode specifically, are primarily 
compressive oscillations. Thus, the molecular cloud turbulence 
can only weakly couple to these oscillations, either through a 
small failure of the incompressible condition or through non- 
linear mode coupling. 

For these reasons we will not consider monopole pulsations 
further, restricting our attention henceforth to oscillations with 
/ > 1. 

2.3.2. Non-radial Oscillations 

Higher-order pulsations are nearly always stabilizing, i.e., 
they increase the allowed surface pressure at which stable 
isothermal cores exist, or alternatively increase the mass that 
a core can support. The exception being the n = 1, / = 1 os- 
cillation which produces a small negative pressure near Ccnfl 
From Equation and the asymptotic behaviors of the ui„i m 
for the non-radial oscillations, it is straightforward to estimate 
the contribution to the surface pressure from the oscillations: 





p 




p 

^osc 




E b 


V ~ 





Po, 



(5) 



where Po is the surface pressure of the unperturbed configura- 
tion. The precise constant of proportionality depends upon ( 
and weakly on the energy spectrum. 

Because for / > the form of u)„i m (R) is well approximated 
by either R' 1 (n> 1) or R~ 3 ' 2 (n = 0), in practice it is not nec- 
essary to fully specify the E„i m . Rather, it is sufficient to define 

E L = E„=o,/>o, m £ "'»> and E osc = Hn>u>a,m E nim, correspond- 
ing to the energy in fundamental oscillations and the p-modes, 
respectively. The resulting excess pressure due to the pulsa- 
tions is then 

E° 



p>o 

J osc . osc 



*osc - 2V ' 3V ' (6) 

where V = 4irR 3 /3 and the difference between the denomi- 
nators arises from the different asymptotic behaviors of the 
uj„i m (R). Note that this is confined to the somewhat narrow 
range from E osc /3V < P osc < E osc /2V. While we will indeed 
specify the full energy spectrum, the distinction between differ- 
ent energy spectra will primarily be in the relative importance 
of£« sc and£>°. 

The particular form of the energy spectrum reflects both the 
mechanism by which the pulsations are excited and their subse- 
quent evolution within the core. Here we consider three differ- 
ent forms for the turbulent energy spectrum, seeking to bracket 
the possible consequences. A natural place to begin is to as- 
sume a Kolmogorov spectrum for the oscillations, i.e., 



-11/3 



if/>0 
otherwise 



where 



k nlm = - v /l+n 2 T/(7+l), 
K 



(7) 



(8) 



so that Ek = 4-nkf jlm E^ lm oc k~ 5 / 3 . Despite the heavy weight upon 
long-wavelengths, this is dominated by n > modes, primar- 
ily due to the dipole and quadrupole oscillations, such that we 
have E^/E® sc ~ 3.5. For concreteness, in FigureQ]the result- 
ing surface pressure is shown by the red line for £ osc = 0.3£b, 
though this scales linearly with E osc /E D . 

4 The fundamental dipole mode, n = & I = 1, does not exist as a conse- 
quence of momentum conservation. 



While not motivated by any particular turbulence model, 
a flat energy spectrum does have the virtue that all realistic 
choices for the turbulent spectra lie between it and the Kol- 
mogorov distributions. Therefore, we also consider, 



nlm 





if«<5and0</<5 
otherwise, 



(9) 



which is flat for n,Z < 5 and vanishing otherwise. The trun- 
cation is arbitrary, here set by numerical convenience. How- 
ever, as argued above, the particulars of the cut-off are not 
important in practice since it is primarily the relative values 
of Eg SC and E>® that matter. As with the Kolmogorov spec- 
trum the latter dominates with E^./E„ sc ~ 3.9 implying that 
the resulting perturbation to the surface pressure will be sim- 
ilar. Indeed, despite the considerable difference in their par- 
ticulars, since E>®/E® sc are similar for both spectra the result- 
ing surface pressures are nearly indistinguishable. The surface 
pressure for the Flat energy spectrum is shown in Figure Q] 
by the blue line (when £ osc = 0.3£b). For both cases we find 
Posc/Po — 0.5 \E osc /Eb\, implying that the fractional change in 
the maximum stable mass of an oscillating isothermal core is 
approximately 0.25 \E 0SC /Eb\. 

Finally, we consider a fundamental-dominated energy spec- 
trum, i.e., one in which only the n = modes are excited, 
with an otherwise Kolmogorov spectrum, as described for E^ lm . 
That is, 



£-11/3 

K nlm 





ifn = 0and/ >0 
otherwise 



(10) 



Aside from presenting an extreme alternative, such an energy 
spectrum will be motivated in Section @] where we discuss the 
evolution of the energy spectrum during core formation. In this 
case £ osc = Eg SC explicitly, and the resulting surface pressure, 
shown by the green line in Figure [T] (again for £ osc = 0.3£bX 
is substantially larger than those produced by the previous 
two spectra. This is a result of both the dependence of P osc 
upon the division of energy between E® sc and as well 
as the departure from the asymptotic expressions of the w„; m 
near £ cr i t for « > 0. As a consequence, in this case we find 
Posc/Po — L0|£ O sc/£b|, corresponding to a fractional increase 
in the maximum stable mass of approximately 0.5 \E osc /E\,\. 

In summary, we have found that the non-radial oscillations 
generally produce an increase in the surface pressure of the 
isothermal cores, with magnitude linearly related to the total 
energy in the pulsations and bounded near by 



0.5 



p 

^ 1 osc 
Pail 



< 10 



(11) 



for an extraordinarily broad class of energy spectra. This cor- 
responds to an increase in the maximum mass of stable cores 
of approximately 



0.25 



E 



< — <0.5 

M crit 



Eb 



(12) 



3. THE STRUCTURE OF MODE-SUPPORTED 
EQUILIBRIUM CONFIGURATIONS 

Thus far we have characterized the effects of oscillations in 
terms of the unperturbed core parameters. However, large- 
amplitude pulsations can dramatically alter the observed ap- 
pearance of a core, with important con sequences for the in - 
terpretation of starless core statistics. In Ke to & Field! (120051) . 
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FIG. 3. — Column density profiles and representative maps for pulsating cores with £ = 14, E osc = 03Eb and a Kolmogorov energy spectrum (Equation|7J. Upper 
Left: Probability distribution of the azimuthally-averaged column density distribution, normalized by the central column density of the critical isothermal core of the 
same mass. The gray-scale is logarithmic, with the floor at a probability density of 0.01 . The inset shows the probability distribution of the central column density. 
For reference the central column density of the static configuration is also shown. On the top axis the £ corresponding to a static isothermal core with the same mass 
is shown. Bottom Left: Isothermal core profiles (dashed red lines, from bottom to top, f = 6, 10, 18, 22, 26 & 30; and solid red line, f = 14.04) and illustrative core 
profiles (green dotted lines) for cores with the central column densities at which the cumulative probability equals 5%, 35%, 65% & 95%, shown explicitly in the 
inset. Center & Right: Example column density maps, corresponding to the profiles shown in the lower left panel. The gray-scale ranges from £/E cr it = to 1.5, 
and subsequent contours are separated in steps of 0.2, with black contours for values less than unity and white above. 



iKeto et al.l (120061) . and lBroderick et al.l (120071) we showed how 
pulsations alter the molecular line profiles and widths. Here 
we explore how pulsations modify the apparent core column 
density profiles, T,(R). In this section we study the range of 
perturbed H(R), finding that generally these are substantially 
asymmetric, frequently devolving into many peaked struc- 
tures, and are biased towards azimuthally-averaged profiles 
that are best fit with super-critical static isothermal configu- 
rations. That is, stable, pulsating starless cores will typically 
appear to be strongly asymmetric and super-critical if fit with a 
static BE model. 



3.1. Generating Column Densities of Pulsating Cores 

To ascertain the typical observational consequences of the 
oscillating cores, we generate an ensemble of pulsating cores. 
We do this for each of the energy spectra described in Section 
I2.3l bv setting the density to 



p(r) = P0 (r) + 5R 




"" Pni(r)Yi„,(r) 



, (B) 



where the po(r) is the unperturbed density, p„i(r)Yi,„(r) are the 
density perturbation eigenfunctions and the $„/,„ are randomly 
determined phases, corresponding to a random realization of 
the particular oscillations. Any realization for which p(r) van- 
ishes anywhere is discarded, though this has little effect upon 
the statistical quantities since the distribution of the various 
observables is nearly identical for the accepted and discarded 
cores. 

We do not vary the E„i,„, nor do we consider different under- 
lying equilibrium configurations, choosing for concreteness a 
stable configuration, near critical, with £ = 14. Since we vary 
only the phases, $„/,„, we are not producing a complete en- 
semble; such an ensemble would necessarily have variations in 
the amplitudes of the oscillations, the total oscillation energy 
and the core masses, in addition to the mode phases. Rather, 
in the absence of strong constraints upon those quantities, we 
restrict ourselves to highlighting the consequences of the pul- 
sations upon the measured densities with fixed amplitudes and 
core masses. 

Once we have generated a p(r) we construct the column den- 
sity map, S(R), by integrating along the z-axis. These are nor- 
malized to the central column density of the static, critical con- 
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FIG. 4. — Column density profiles and representative maps for pulsating cores with £ 
see Figure[5] 



: 14, £ osc = 0.3£b and a Flat energy spectrum (Equation|9j- For further details 



figuration, E cr ; t . While the center of the unperturbed core is 
always located at the origin, we define the "central" column 
density, E c , to be max [E(R)], corresponding to how the center 
of the core would be defined in practice. Thus, configurations 
which have £ c /£ cr it > 1 appear to be super-critical. 

Finally, we produce azimuthally-averaged column density 
profiles, Ti(R), by averaging in concentric annuli about the 
peak column density. Since the perturbations can both raise 
and lower E(R), this observationally motivated definition of E c 
necessarily biases us towards higher column densities, and thus 
towards identifying pulsating cores as closer to super-critical 
than the are. 

We repeat this process many times, building up a statistical 
description of the pulsating core column densities. The results 
are summarized for E osc = Q.3E D in Figures [3}|5] and are dis- 
cussed in detail below. 

3.2. Central Column Densities 

The probability of measuring a given central column den- 
sity is shown in the inset of the upper-left panel of Figures 
|3]-0 In all cases this probability is broadly distributed about 
E c /E cr it = 1. As with their dynamical consequences, the distri- 
butions arising from pulsations with Kolmogorov {E^ lm ) and 
Flat (£,f/,„) energy spectra are quite similar. Both probabil- 
ity distributions are centrally peaked at a value significantly 
larger than unity, E e /E cr i t ~ 1.2, corresponding to a £ ~ 21, 



with only a roughly 5% probability of finding £ < £ cr j t (seen in 
the cumulative probability distribution in the inset of the lower- 
right panels of Figures |3f|5j. Furthermore, these can appear to 
be extraordinarily super-critical, exhibiting £'s of up to 28 for 
E osc = 0.3£b- 

Differing from these qualitatively, the distribution arising 
from the fundamental-Kolmogorov energy spectrum (E^) is 
double peaked. This is due to the dominance of the fundamen- 
tal quadrupole mode (the dipole has no fundamental). It is also 
nearly symmetrically distributed about E c /£ cl -i t = 1, somewhat 
more narrowly than with the Kolmogorov and Flat spectra, and 
with < Ccrit) — 50%. Despite this, the double peaked na- 
ture of the distribution implies that most cores will appear to 
have £'s significantly larger or smaller than £ cl i t , reaching up to 
roughly 17 for E osc = 0.3E D . 

3.3. Column Density Profiles 

The results from the central column densities are borne out 
by the azimuthally averaged column density profiles. As be- 
fore, the Kolmogorov and Flat energy spectra produce quan- 
titatively similar results. The probability distributions of 
£(/?) /E cr i t are shown in the leftmost two panels in Figures[3}|5] 
Over-plotted in red within the bottom left panels are a number 
of static isothermal core profiles, for a number of £. £(/?) for 
the critical case is shown by the solid line, with profiles mov- 
ing upward and downward in £ increments of 4. In particular, 
the £ = 18 and 22 profiles do a fair job of fitting the averaged 
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FI G. 5. — Column density profiles and representative maps for pulsating cores with f = 14, £ osc = ().3£b and a fundamental-Kolmogorov energy spectrum (Equation 
110) . Note the change in scale in the left-panel insets. For further details see Figure[3] 



column density profiles, while those with significantly smaller 
( do not. For reference, profiles associated with randomly cho- 
sen examples with central densities at the 5%, 35%, 65% and 
95% cumulative probabilities are shown by the green dotted 
lines (in each Figure, these are associated with the maps shown 
on the right). 

The fundamental-Kolmogorov energy spectrum produces 
similar column density profiles. As anticipated by the £ e /£ CI j t 
distribution, these are centered upon the critical profile, with 
both low and high-£ profiles represented. In this case the av- 
eraged profiles appear to fit the static-core profiles more accu- 
rately, though this may be due to the decreased deviation from 
the critical core profile. 

This suggests that if observations of pulsating cores are com- 
pared to a static BE model, the observed cores will quite com- 
monly appear strongly super-critical with f > ( crit even though 
the cores are in no imminent danger of collapse. 

3.4. Column Density Maps 

Despite the approximate equivalence between the average 
column density profiles of stable pulsating cores and static 
super-critical cores, the full two-dimensional maps of pulsat- 
ing cores can be dramatically asymmetric. Example column 
density maps, randomly chosen, are shown in the center and 
right columns of Figures [3]-[5] corresponding to the same con- 
figurations for which the azimuthally averaged distributions are 
shown in the bottom left panel of each Figure. These are repre- 



sentative of the configurations appearing at the 5%, 35%, 65% 
and 95% cumulative probabilities, giving some idea of how 
typical and, perhaps, rare cores appear. 

Again the column density maps associated with the Kol- 
mogorov and Flat energy spectra are similar, though maps as- 
sociated with the former are typically more centrally concen- 
trated. In roughly 20% and 10% of cases, for the Kolmogorov 
and Flat spectra, respectively, the column density maps are 
double-peaked. 

With fewer degrees of freedom, and smaller net den- 
sity variations for the same E osc , it is not surprising that 
the fundamental-Kolmogorov spectrum shows less small-scale 
structure. In most cases the column density maps are only 
weakly non-circular. Driving E osc up to 0.6£b produces 
roughly comparable variations in E f /E cl i t , and correspond- 
ingly larger deviations from cylindrical symmetry, typically 
with aspect ratios larger than 2. 

Nevertheless, for all three energy spectra considered it is pos- 
sible to produce manifestly asymmetric structures, comparable 
to those observed. As a consequence, strong departures from 
symmetry in the column densities are not generally evidence 
for unstable structures. 

4. EVOLUTION OF MODE-SUPPORTED STARLESS CORES 

In addition to the excitation mechanism, the evolution of the 
underlying core properties and the non-linear damping of oscil- 
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lations can imprint themselves upon the pulsation energy spec- 
trum. Ultimately, the relative importance depends upon the re- 
lationship between the timescale associated with each process. 
Here we present some observationally motivated constraints 
and discuss the consequences for the stability and appearance 
of starless cores. In particular we motivate a picture in which 
the oscillations are inherited from the parent molecular cloud 
turbulence, and the fundamental modes are amplified at the ex- 
pense of shorter wave length during a slow contraction phase 
resulting in a fundamental-mode dominated spectrum. Sub- 
sequent non-linear decay of the pulsations then eliminates the 
additional turbulent support leading to collapse and star for- 
mation in those cores that are super-critical with respect to the 
static BE model. 

4.1. Timescales 

The relevant timescales for the oscillating starless cores are 
the evolution timescale of the core, T, the pulsation periods, 
Priim, lifetimes, r„/ m , and excitation timescale t e . Even with- 
out knowing the excitation mechanism we can begin to place 
these into a hierarchy based upon the observation that cores 
with large-amplitude, long-period pulsa t ions exist. In addition 
to t he observations oflLada et alj ( 12003b , iRedman etail ( 120061) 
and lAguti et aTl (12007 ) mention ed in the introduction, the ob- 
servations of LSohn et alJ J200l find that about 2/3 of all their 
85 observed cores show clear evidence of either expansion or 
contraction. Most of the remaining 1/3 show complex spectral 
line profiles that also indicate internal motions but are diffi- 
cult to characterize. The presence of coherent oscillations im- 
mediately implies that P„i, n < T, t„i„ u t e . That is, for coherent 
oscillations to exist the core must persist in a fixed state suffi- 
ciently long for the pulsation to coordinate the large-scale mo- 
tions throughout the core, and the pulsation amplitudes cannot 
change significantly during this process. Violations of either of 
these would result in motions that appear stochastic (microtur- 
bulent) instead of coherent. 

That the observed pulsations have E osc ~ E\, implies r„;„, < t e . 
Were the observed modes the result of a prolonged excitation, 
or many individual excitation events, we would expect with 
near unit probability that a given core will have had a sequence 
of excitations resulting in £o SC ^s> E\,, at which point the core is 
torn apart, in direct conflict with the observation of numerous 
cores with substantial internal motions. In principle, this can be 
prevented by limiting the energy transfer during each excitation 
event. However, in that case the development of large-scale 
pulsations is itself rare, again in direct conflict with observa- 
tions. Thus, as long as the excitation mechanism is uncorre- 
cted with the pulsations themselves, i.e., there is no oscillation 
"feedback", a given mode must have, on average, originated 
from only a single excitation event. 

Finally, we have T < t„/,„. As the pulsations damp away, the 
support arising from the oscillations is lost, forcing the core 
to evolve towards higher central densities and then ultimately 
towards collapse. Thus, for cores in which the pulsations are 
dynamically important the mode lifetimes set an upper limit 
upon the core evolution time; a core may evolve on shorter 
timescales but it will evolve as modes decay. 

Thus, in summary, the very presence of starless cores with 
dynamically significant pulsations gives the following hierar- 
chy among the relevant timescales: 



Kolmogorov 



Flat 



Pnlm <T <T„l m < t e 



(14) 
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FIG. 6. — Evolution of the pulsation energy spectra during adiabatic core 
collapse. The dotted lines show k^i m E„i m /E ox for the Kolmogorov (left) and 
Flat (right) energy spectra, assumed to be inherited at (J = 1 .25, as functions of n 
for various /. The solid lines show the adiabatically evolved spectra at £ = 14, 
with in all cases the fundamental (n = 0) modes gaining substantially more 
energy than the other oscillations. In the case shown, by f = 8 the spectra have 
nearly attained their final state, implying that even for substantially sub-critical 
cores, the pulsation energy spectrum can be dominated by the fundamental 
modes. 

4.2. Energy Spectrum from Inherited Turbulence 

Given that T < t e , we are justified in imagining that pulsa- 
tions on super-critical cores are relics of the core formation 
process. The natural mechanism for the excitation of the oscil- 
lations is then the turbulence within the parent cloud, inherited 
at the formation of the core. This does not mean that the en- 
ergy spectrum of an observed core is always the same as the 
parent cloud. If cores are formed by a prolonged collapse, the 
evolution of the underlying equilibrium configuration modifies 
the energy spectrum of the oscillations. Since P„i„, < T, this 
process is slow in the sense described in Section l2~2l and again 
we may use the adiabatic invariant E„i m /ui n i m to estimate the 
evolution of the turbulent spectrum. If this formation occurs at 
constant temperature and mass, the mode periods are functions 
of ( alone, with the dependencies implied by Figure [2] Thus, 
the evolved energy spectrum is given by 



EnlmiC\ Co) = E, 



K/F ^>nlm(Q 

w„;,„(Co) 



(15) 



This simplifies dramatically the determination of the conse- 
quences of mode decay and core evolution. 



where £o is the center-surface density contrast at the time that 
the molecular cloud turbulence was inherited. The spectrum 
of supersonic turbulence in larger scale molecular clouds is not 
well understood. Observations indicate that it is close to Kol- 
mogorov, in the range of a ~ L" 1 / 3 to a ~ L" 1 / 2 where a is the 
observed velocity dispersion and L the length scale. For con- 
sistency with our previous analysis, we will consider both our 
Kolmogorov and Flat energy spectra as the inherited spectra. 

For the same reason that the fundamental modes (n = 0) 
and p-modes (n > 0) produce different contributions to P osc , 
an adiabatic contraction results in a fundamental-heavy energy 
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spectrum if £o is sufficiently small. This is shown at £ = 14 
in Figure [6] for both the Kolmogorov and Flat energy spectra 
with (o = 1.25. During the collapse phase, E^/E^ c has de- 
creased from 3.5 and 3.9 to 1.3 and 1.5 for the Kolmogorov 
and Flat energy spectra, respectively. Thus, while in neither 
case the fundamental modes dominate, in both the energy in 
the fundamentals and p-modes are comparable. This suggests 
that neither E^ m nor EL may be applicable, but rather some 
combination of the two is appropriate for the observed pul- 
sating starless cores. This also provides a natural explanation 
for the dominance of low-order modes, and in particular low- 
order quadrupole modes, as implicated by sel f-absorbed asym- 
metric molecular line spectra in Barnard 68 (|Lada et aTlfe OCB ; 
iKeto et al.ll2006llRedman eTaIll2006llBroderick et al.ll2007l) . 

4.3. Pulsation Damping and the Onset of Instability 

In Section [2] we showed how pulsations can support super- 
critical isothermal configurations. However, in the absence 
of continual excitation, these oscillations eventually decay. 
The manner and timescale in which large-amplitude pulsations 
damp depends upon the underlying core parameters and en- 
vironment as well as the detailed structure of each individ- 
ual oscillation and their amplitudes. In general the lifetimes 
are different for each mode. In our previous papers we have 
estimated life times for a few particular cases. Figures 2 
and 5 of iBroderick et all d2007l) estimate a decay time of 2 
Myr for the 122 and 1 1 1 modes by mode-mode coupling. In 
IBroderick et alj d2008l) . we estimate decay times of 1 Myr and 
0.4 Myr by transmission to the larger scale ISM for isolated and 
embedded cores respectively. These examples are too few to 
define the decay times precisely; therefore, for simplicity, we 
assume that all pulsations decay exponentially with the same 
lifetime, r ~ 10 6 yr, chosen to be consistent with our previous 
numerical simulations of the non-linear mode-mode and core- 
environment couplings (IBroderick et alJl2"007ll2008h . 

The consequences of mode damping are shown explicitly in 
Figure [7] Here P(R) curves are shown for a number of mode 
amplitudes, where the energy spectrum is that resulti ng fr om a 
slow initial contraction period, described in Section 14.21 Ini- 
tially, these are normalized such that 50% of the binding en- 
ergy is in oscillations at £ = 14 (dark-blue line), and then at 
f/r = 0.5, 1.0, and 1.5. Since the exterior pressure is constant, 
we show a number of points with a fixed, super-critical surface 
pressure, illustrating the onset of instability during Pulsation- 
decay driven evolution. The associated evolution in which in 
this case is simply proportional to the central density, is shown 
in the inset, with the colored points corresponding to the curves 
in the main Figure. 

Prior to the times shown, we assume the cores were forming 
with T < t, presumably during a phase of prolonged contrac- 
tionQ At some point, due to both the increase in the dynamical 
importance of the pulsations and the decrease in r due to the 
increased mode amplitudes and central density, the oscillations 
begin to dominate the formation processes, and the cores en- 
ter into a pulsation-driven phase. During this time cores can 
derive a substantial portion of their pressure support from the 
oscillations themselves, and their damping drives changes in 
the underlying configurations on timescales comparable to r. 
Finally, for super-critical cores, once the energy remaining in 

5 Note that r is also a function of the core density, and nonlinearly dependent 
upon the dynamical timescale of the core. Thus, there always exist R for which 
the mode-damping timescales are larger than the free-fall timescales. This is 
not the case for the mode periods, however. 
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FIG. 7. — Surface pressure as a function of radius for an pulsating isother- 
mal gas sphere as the oscillations decay. Initially, at £ = 1.25, the gas has a 
Kolmogorov energy spectrum, evolving as described in Section l4~2l thereafter, 
and normalized such that at t = iJosc = 0.5£b at £ = 14. The pulsations energy 
decays exponentially in time, with decay constant r. The dark-blue, light-blue, 
light-red and red lines show the resulting surface pressure at t = 0, 0.5r, l.Or 
and 1.5t, respectively. The black points show the equilibrium configuration at 
a fixed, super-critical surface pressure. The inset shows time evolution of the 
central density. For reference, the points corresponding to the curves shown in 
the main figure are also plotted. Finally, the relative location of the Formation 
and Collapse phases are shown. 

pulsations has reached some critical value, oscillations can no 
longer support the core, leading to unstable collapse. Thus, we 
may think of the evolution of a pulsating starless core as con- 
sisting of at least three phases: Formation, Pulsation-supported 
and Collapse. 

4.4. Implications for the Core Mass Function 

We have already used the existence of large-amplitude oscil- 
lations to empirically determine the relationship between the 
mode excitation timescale and the other relevant timescales. 
However, these timescales are not generally independent of the 
core mass, and may be violated for cores more massive than 
those observed, with consequences for the core mass function. 
Here we consider these within the turbulent excitation picture. 

The starless cores, or their progenitors, are constantly being 
buffeted by the turbulence within the parent molecular cloud. 
The ability of a core to survive long enough to potentially form 
a star requires that the turbulence not excite pulsations with 
sufficient amplitudes to tear the core apart, either by a single 
catastrophic excitation or a sequence of smaller excitations. 
We begin by noting that in the absence of oscillations there 
is a well defined M(R), shown in Figure [T] and at low masses 
roughly oc,R~ 3 / 2 . When dynamically significant oscillations are 
present, the core mass is a function of radius and E osc . Never- 
theless, that low mass cores are generally large and vice versa. 
As a consequence, the binding energy of cores, comparable 
to E h (R) ~ GM 2 /3R (within 30% at all masses), is a strong 
function of core mass, and in particular considerably larger for 
massive, compact cores. 
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In contrast, the energy in turbulent motions (for Kolmogorov 
turbulence) on the core scale is E t oc R 2 ^, which decreases 
with size. This implies that for some radius, (or mass, 
Mo), which depends upon the strength of the turbulence, we 
have Eb(Ru) = E t (Ru). Since the turbulent motions fluctuate on 
timescales comparable to the circulation time of a single eddy 
(i.e., 2ir times the sound crossing time of the core), which is 
itself comparable to the core dynamical timescale, we expect 
to see at least one such excitation event for every core. For 
large, low-mass cores we have E\, < E t , and we expect the core 
to be torn apart when this occurs, placing a lower-limit upon 
the mass of long-lived cores. In contrast, for compact, massive 
cores we have £\, > Eu an d external turbulence is no longer 
able to excite large-amplitude motions within the core leading 
aga in to the inherited turbulence picture described in Section 

M 

The fact that E^jE^ ~ (R/Ro)~ 14 / 3 is such a strong function 
of R implies that the range in core radii about Rq, and thus mass 
about Mo, in which large-amplitude pulsations can be excited 
is narrow. In particular, 
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Cores much less massive than Mo are torn apart, while those 
much more massive than Mo never develop the observed large- 
amplitude pulsations. The latter issue is moderated somewhat 
if an adiabatic contraction phase occurs during the formation 
process, i.e., if the pulsation amplitude are representative of 
those generated when the core was considerably larger. Never- 
theless, the former places a hard constraint upon the minimum 
mass of starless cores. 

5. CONCLUSIONS 

The onset of instability in isothermal gas spheres depends 
upon the dynamical state of the object. While this has long 
been known, usually within the context of effective turbulent 
support, how this occurs in practice depends upon the oscil- 
lation mode energy spectrum, and therefore the mechanism 
by which pulsations are excited. In principle, large-amplitude 
breathing modes can induce instability for masses well below 
C > 4 the static stability limit, a fact that is reflective of the na- 
ture of the instability. However, in practice, doing so requires 
extraordinarily non-linear mode amplitudes, at which point the 
dynamical effect is unclear. In contrast, higher-order pulsa- 
tions are generically supportive, increasing the maximum sta- 
ble mass significantly. For oscillation energies comparable to 
those required to explain the self-absorbed, asymmetric molec- 
ular line profiles typically observed in starless cores, the maxi- 
mum stable mass is about 10-30% higher than the static case. 

If the modes are excited during the core's formation, pre- 
sumably inherited from the turbulence of the parent molecu- 
lar cloud, followed by an adiabatic contraction, it is natural 
to expect low-n oscillations to dominate the energy spectrum. 
That is, even prior to mode decay via mode-mode or mode- 



environment coupling, energy from the collapse is preferen- 
tially shunted into the long-wavelength pulsations. Thus, it is 
not unreasonable to expect starless cores to pulsate predomi- 
nantly in low-order multipole modes, with / > 2. 

The presence of large-amplitude oscillations can dramati- 
cally alter the morphology of the column density maps of star- 
less cores, mimicking the unstable configurations of static BE 
spheres. For turbulent energy spectra, dynamically significant 
pulsations dramatically bias the maximum column densities to- 
wards larger values, with the average exceeding that associated 
with the critical configuration of static BE models by factors of 
order unity. This is born out by the azimuthally-averaged col- 
umn density profiles, which are very similar to isothermal gas 
profiles associated with £ above the critical value. Despite this, 
the two-dimensional column density maps are strongly asym- 
metric, featuring in some cases multiple peaks, similar to those 
observed in practice. 

The ultimate decay of the oscillations, facilitated by non- 
linear mode coupling or by coupling to the surrounding molec- 
ular gas, results in the collapse of super-critical cores. For slow 
decay (lifetimes exceeding the mode periods) this results in a 
sequence of quasi-stable isothermal gas spheres, and the evolu- 
tion may be determined explicitly. This assumes that the oscil- 
lations are not continually re-excited by small scale processes 
within or near the starless cores themselves. However, the very 
existence of a significant number of pulsating cores, with mode 
energies being a considerable fraction of the cores' binding en- 
ergy, argues against this. 

Taken together, these suggest a picture of star formation, 
framed entirely within a hydrodynamical paradigm. Self- 
gravitating cores condense out of the turbulent sea within 
molecular clouds, born with large-amplitude pulsations. These 
pulsations drastically perturb the core shape, mimicking non- 
equilibrium structures as well as initially supporting the cores 
against collapse. As the oscillations decay the core becomes 
unstable resulting in star formation. In this scheme the rate of 
star formation is set not by ambipolar diffusion of a large scale 
magnetic field, but by the non-linear decay of the core pulsa- 
tions. 

We should emphasize that while we have framed this new 
picture entirely in terms of hydrodynamics, small scale tangled 
magnetic fields in equipartition with the turbulent energy are 
still allowed and might be required. The model of oscillations 
as perturbations does not allow amplitudes large enough to pro- 
duce the supersonic motions occasionally seen in a few cores. 
These are the exception rather than the rule and might not be 
pulsations at all, but rather inward motions as the core center 
transitions to free-fall collapse. However, if these motions are 
sub-Alfvenic pulsations, they may be subsumed into an equiva- 
lent magnetohydrodynamic formulation, in which we consider 
only small-scale magnetic fields. In this case, it would be the 
decay of the magnetohydrodynamic pulsations, not ambipo- 
lar diffusion, that sets the star formation rate. Nevertheless, 
as empirical evidence accumulates for the existence of large- 
amplitude oscillations in starless cores, it is clear that their dy- 
namical effects upon the cores themselves cannot be ignored. 
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